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SUMMARY 

An improved finite difference procedure for determining the natural fre- 
quencies and mode shapes of tapered cantilever beams undergoing uncoupled vi- 
brations is presented. Boundary conditions are derived in the form of simple 
recursive relations involving the second order central differences. Any ap- 
proximation error resulting from this process is discussed. Results obtained 
by using the conventional first order central differences and the present 
second order central differences are compared, and it is observed that the 
present second order scheme is more efficient than the conventional approach. 
An important advantage offered by the present approach is that the results 
converge to exact values rapidly, and thus the extrapolation of the results is 
not necessary. Consequently, the basic handicap with the classical finite 
difference method of solution that requires the Richardson's extrapolation 
procedure is eliminated. Furthermore, for the cases considered herein, the 
present approach produces consistent lower bound solutions. 


INTRODUCTION 

Several methods of solution of beam vibration problems, which can be 
broadly classified as either belonging to the continuum model approach or to 
the discrete model approach, have been published. Principles of the minimum 
potential energy, or complimentary energy, the Reissner mixed method and the 
Dean and Plass principle have been used in the continuum model approach. It 
has been established that the potential and complimentary energy principles 
give an upper bound to the solutions while the Reissner method and Dean and 
Plass method may produce upper bound solutions depending on the type of formu- 
lation and choice of shape functions. However, bounds for these mixed methods 
have not been established theoretically. Each of these methods has its inher- 
ent advantages (refs. 1 to 3). In the discrete model approach, the Holzer, 
Stodola, polynomial frequency equation, Myklestad, and finite element methods 
are all well known. Solution of the equations of motion in the continuum 
model approach is possible under certain conditions by the application of the 
Galerkin process, the collocation method, and the finite difference method. 
Among these techniques, the Galerkin process is known to be equivalent to the 
Rayleigh-Ritz methods and generates upper bound solutions identical to those 
obtained by the Ritz method, provided that similar shape functions are used in 
both techniques. It has been observed that the accuracy of the results ob- 
tained by using the collocation method depend on the choice of the collocation 
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points and their location (ref. 4). The finite difference method has attrac- 
ted considerable attention and several contributions exist which make use of 
the first-order forward, backward, or central difference schemes or their com- 
binations (refs. 5 to 12). In almost all the works, it has been observed that 
the finite difference method is subject to relatively slow convergence with 
mesh refinement, although Richardson's extrapolation procedure (ref. 13), when 
applied to two or three successive iterations with different mesh sizes, can 
produce results which may be close to the exact solutions. However, such an 
extrapolation procedure requires that the convergence of the results be mono- 
tonic, and the extrapolated result may not necessarily give a bound. 

Relatively few works exist which deal with the application of higher or- 
der finite differences. Greenwood (ref. 14) used first order and second order 
finite difference schemes, which produce truncation errors of the order O(h^) 
and 0(h'+), respectively, for the analysis of uniform beams in flexure and for 
another case having uniform breadth but varying depth with fixed-free end con- 
ditions. The fourth order differential equation for the uncoupled vibration 
was transformed into four first order equations, the slopes and shearing for- 
ces were evaluated at half integer stations while the deflections and bending 
moments were evaluated at integer stations. Central difference approximation 
of order O(h^) with staggered stations was used. One-sided approximations 
were used to satisfy the boundary conditions wherein all the stations encoun- 
tered were inside the beam and fictitious stations outside the beam were not 
required. As an alternative approach, another complicated, but symmetric, 
method of representing the boundary conditions was also illustrated. It was 
concluded that one-sided approximation gave better results while the symmetry 
assumptions gave results of moderate accuracy. It is interesting to note from 
Greenwood's results that, for the tapered beam case, the O(h^) approximation 
yields better accuracy than the higher order O(h^) approximation. It was 
stated in reference 14 that the lumping of nonuniform mass caused the largest 
loss of accuracy in the O(h^) approximation. 

Gawain and Ball (ref. 15) presented finite difference formulae having 
consistent errors of O(h^) by representing the function with a truncated 
power series at the boundary in such a manner that the boundary conditions are 
satisfied. A consistent error of O(h^) was obtained for the error estimate. 
However, higher order central difference expressions having truncation errors 
of order O(h^) have not been reported. 

In what follows, a second order central difference approach is presented 
which eliminates most of the problems discussed above. Clamped free beams are 
analyzed for axial, torsional, and flexural vibrations. Boundary conditions 
are enforced which eliminate the fictitious stations outside the beam. These 
conditions are obtained from the simple and logical extensions of the first 
order theory. Any approximation error encountered in this process is dis- 
cussed. It will be shown that the present theory produces accurate results 
with rapid convergence and, consequently, the extrapolation procedures needed 
in the classical first order central difference theory can be successfully ob- 
viated with the present improved formulations. 

Professor A. W. Leissa provided facilities at Ohio State University from 
March to April 1983, during which time the uniform beam cases were solved. 

The results of the uniform beam cases are reproduced in this report for 
completeness. 
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Subscripts: 


area at any section 
breadth of beam 
torsional rigidity 

derivative with respect to z 

second derivative with respect to n 

Young's modulus 

modulus of rigidity 

length of each elemental beam segment 

second moment of area about flexible plane 

polar moment of inertia about centroid 

length of beam 

number of beam segments 

natural radian frequency 

time 

depth of beam 

dynamic displacement in longitudinal direction 

dynamic displacement in flexible direction 

coordinate measured along longitudinal direction of beam 

breadth taper parameter 

pretwist over length of beam 

depth taper parameter 

axial fractional length 

dynamic torsional displacement 

mass density 


-1, -2, n+1, 
n+2, n+3 
Superscripts: 


arbitrary station 

fictitious stations outside beam domain 


differentiation with respect to z or n 
differentiation with respect to time 
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ANALYSIS 


Axial Vibrations of Tapered Cantilever Beams 

The governing differential equation for free axial vibrations of a canti- 
lever beam is 


. d 


to 


(eA ♦ PP^W - 0 (1) 

For a beam with breadth taper 3 and depth taper 6, equation (1) reduces 

aw" + bw' + cp^ = 0 (2) 


where 


a = (1 - 3n)(l - 6n)E/L‘ 


b = -E[B(1 - «n) + 6(1 - Bn)]/L‘ 


c = p(l - 3n)(l - 6n) 


n = z/L, 0 < n < 1 




(3) 




The boundary conditions for the cantilever beam reduce to 

w(0) = w' (1) = 0 


(4) 


Solution by first order central differences. - Substituting the central 
differences for the derivatives of w given in the appendix in equation (2) 
and simplifying, one can write the following equation for any arbitrary station 
i of the beam as follows: 


^i'^i-1 ^i'^i ^i^i+1 = i = 0, 1, . . n 


(5) 


In the preceding equation 


A^. = a. - b.h/2; B. = -2a.; C. = a. + b.h/2; D- = -c.h^; 
a. = (1 - 3 ih)(l - 6ih)E/L^; b. = -E[b( 1 - 6ih) + 6(1 - Bih)]/L^;[" (6) 

c^. = p(l - Bih)(l - 6ih) 
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The boundary conditions in terms of the first order central differences are 

“o = 0: “n - ® <”■ '•nn - Vl 

where n represents the station located at the free end of the cantilever. 
Equation (5) together with the boundary conditions given by equation (7) leads 
to the frequency equation In the form of n-equations for i = 1, 2, . . n: 
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Solution by second order central differences. - The boundary conditions 
for the fixed-free case of a beam in axial vibration given by equation (4) in 
terms of second order central differences (see appendix) are 


"o = "n ' -Rh <"n-2 - ®"n-l * 


n+1 


- V2^ " ^ 


(9) 


Assumption of the symmetry condition (ref. 14) at the free end for the func- 
tions (wn_i, Wn+i) and (wn_ 2 » Wn+ 2 ) leads to 

“n-l " "nn; ^-2 ' «n+2 ' 


and the fictitious stations w^+i and Wp +2 can thus be eliminated. As can 
be seen from equation (7), the physically consistent condition, Wf^_j = as 

given by the first order theory is extended to cover the fictitious station 

Wn +2 also. Another important observation can be made by inspecting equa- 
tion (5). If this equation is evaluated for the station i = 0, the built in 

end, one obtains the condition 




B.W- + C-W, = 
00 0 1 


^ 00 


( 11 ) 
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which implies that 



1 - 



(e + «)h 

2 

(e + 6)h 
2 


Wi 


( 12 ) 


yielding w_i = wj for a uniform beam. Substituting the second order fin- 
ite difference expressions for the derivatives given in the appendix in equa- 
tion (5) and simplifying, one obtains the following equation for any arbitrary 
station i: 


* CjW, ♦ 


“i"in * 


- PlP^x,- 


(13) 


In the previous equation 


A. = (b.h - a.)/12h^; B- = (16a. - 8b.h)/12h^; 
C. = -30a^/12h^; D^- = (16a. + 12b.h)/12h^; 




> ( 14 ) 


= -(a^. + b^.h)/12h^; = -p(l - 6ih)(l - Bih); 

a. = E(1 - Bih)(l - 6ih)/L^; b. = -E[b(1 - fiih) + 6(1 - Bih)]/L^ 


If equation (13) is evaluated at n = 0, it can be shown that 

w_i = -(Dq/B^)wj = - |[1 - h(B + 6)/2]/[1 + h(B + 6)/2]| Wj 

''-2 = -(^o/^^‘'2 = - (C1 - f’(e «)V[1 + h(B + 6)]| w^ 

and it can be seen that the first of equations (15) is identical to the corre- 
sponding expression given by the first order central differences. 

Making use of equations (10) and (15) and evaluating equation (13) at each 
station i = 1, 2, . . ., n results in the following frequency equation: 
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= p^|(F^w^)(F 2W2) . . . (F^w^)} (16) 


Torsional Vibrations 

Torsional vibrations of tapered cantilever beams: second order finite 

differences. - The governing equation for free torsional vibrations of a beam 
is given by 

The boundary conditions for a fixed-free beam are 

e(0) = 0; e '(L) =0 (18) 

For a tapered cantilever beam of rectangular cross section with width taper 6 
and depth taper 6, torsional rigidity C, and polar moment of inertia Ip 
take the form 


7 


(19) 


C =— 


ez)(l - 62)^ - 


TT 0 




N=l,3,5 L 


1 N, •>o(l - 


T t^u - 6Z) 


T 0 0 

=~TT 


(1 - Bz)(l - «z)[bQ(l - ez)^ + t^(l - 62 )^] 


Substituting equations (19) in equation (17), performing the necessary differ- 
entiation and writing in nondimensional form, one obtains the following equa- 
tion in terms of second order finite differences for any station i: 


+ E.9.^2 = F.p e. 


( 20 ) 


and 


A. = S1(S2 h - S3); B^. = Sl(16 S3 - 8h S2); C. = -30 SI S3; 
D. = Sl(16 S3 + 8h S2); E^. = -S1(S3 + h S2); F^. = pi 


SI = l/12h' 


S2 = ^^ir-36(l - Bih)(l - 6ih)^ - p(l - 6ih)^] 


3L^ L 
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N=1,3,5L 


-S4' 

e - e 

m5 I «S4 . ^-S4 
N \ e + e 
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2iTb (6 - b) 


t^(l - 6ih)2N"^ 


(e^' - e-S') 


S3 = 


Gb t^ 
0 0 

3L^ 


t- 


Bih)(l - 6ih)^ - (1 - 6ih)^ 

IT 
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Ln^ (e^4 , ^-S4j 


S4 = 


Nnb^(l - Bih) 
2tQ(l - 6ih) 
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By proceeding on the same lines as presented for the case of axial vibration 
and noting that the boundary conditions in terms of the finite differences can 
easily be written analogously to those given in equations (10) and (15), one 
can develop the frequency equation which will be identical to the one presented 
in equation (16), with A-j, B-j, . . ., F-j defined by equation (21). 

Torsional vibrations of pretwisted-tapered cantilever beams. - When a 
beam is pretwisted, an increase in the torsional rigidity takes place due to 
the inclination of the blade fibers in addition to the fiber bending effects. 
Carnegie (ref. 16) derived a correction factor for fiber bending and pretwist 
effects for thin rectangular cross section blades. So far, no rigorous mathe- 
matical solutions are available for the fiber bending effects of tapered 
blades. Among the recent contributions on the torsional vibrations of pre- 
twisted blades, at least two works require careful consideration. Duggan and 
Slyper (ref. 10) observed that the boundary conditions adopted by Carnegie 
(ref. 16) do not yield accurate solutions for low aspect ratio blades. They 
drew attention to the boundary condition discussion by Barr (ref. 16) and used 
a modified set of boundary conditions to obtain the torsional frequencies for 
low aspect ratio blades. Kaza and Kielb (private communication) reported on 
the case of rotating pretwisted cantilever blades, allowing for the effects of 
warping rigidity. The torsional equation of motion was derived, and the bound- 
ary conditions were established through a variational formulation (Hamilton's 
principle). They studied the effect of structural warping and of inertial 
warping for blades having wide ranges of aspect ratios and rotational speeds. 
Results were generated by using the Galerkin method. Closed form exact solu- 
tions were obtained for the nonrotating case. These equations were used by 
Subrahmanyam and Kaza (unpublished data) with a first order finite difference 
method of solution for uniform pretwisted blade cases of small and large aspect 
ratios; close agreement between the two approaches was shown. 

If warping is included, the torsional equation of motion will be of fourth 
order with variable coefficients, and incorporation of the effects of taper 
will be more involved. On the other hand, for thin rectangular blades of large 
thickness and aspect ratios, Carnegie's formulation, which leads to a second 
order equation of motion and to an appropriate set of boundary conditions, will 
be adequate, and the equations developed in the preceding section for untwisted 
beams in torsion can be used with only slight modifications. 

If the additional torsional rigidity due to pretwist is incorporated, the 
net torsional rigidity (neglecting structural warping effects) can be written 
as 

C^ = C + C^ (22) 

where C is given by equation (19) and the increase in torsional rigidity due 
to pretwist, over and above that of St. Venant, is 
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(23) 




Et b" 


5 = yZ/L = yn 




The frequency equation for the case of a pretwisted tapered blade will be ob- 
tained by replacing S2 and S3 by S2 and S3, which are defined as 

Et b^ ^ 

S2 = S2 - B2)^ + 56(1 - 6z)(l - 6z)^l 

ISOL'^ L -* 


_ Et b^ 5 j- 

S3 = S3 + — ^ y‘^(l - 6Z)^(1 - 6Z) 
180L^ 


(24) 


Flexural Vibrations of Tapered Cantilever Beams 

The differential equation for flexural vibration of beams neglecting shear 
deflection and rotary inertia effects is of the form 


2 

^.EI - App^y = 0 

dz"^ V dz*^ 


(25) 


which can be written in the following nondimensional form 

^ (eI - pAL^p^y = 0 

dn V dn 


(26) 


The boundary conditions for a fixed-free beam can be reduced to 


y = 0 and = 0 at n = 0 




.2 .3 

— ^ = 0 and ■ ' A =0 at n = 1 
dn dn ^ 


(27) 


Equation (26) is rewritten as fallows for a tapered beam of rectangular cross 
section: 
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( 28 ) 


a^+b-^+c-^ = dp^y 
dn dn dn 


where 


a = (1 - 6n)(l - 5n)^ 
b = -2 b( 1 - 6n)^ - 65(1 - Bn)(l - 6n)^ 
c = 665(1 - 5n)^ 65^(1 - Bti)(1 - 5n) 

d = 12pL^(l - 6n)(l - <5^)/Et^ 






(29) 


Solution by first order central differences. - The boundary conditions 
represented by equation (27) when written in terms of first order central dif- 
ferences give 


^ = 0’ ^-1 = yp ^n+1 = ^^n - Vl I 

^n+2 = ^n-2 " ^^n-1 ^ ^^n j 


(30) 


By using the finite difference expressions for the derivatives from the appen- 
dix, incorporating the relations given by equation (30) and following on the 
lines described earlier, the frequency equation for the flexural vibrations 
can be developed easily and is identical to that given by Carnegie and Thomas 
(ref. 11). 

Solution by second order central differences. - Substituting the second 
order finite difference equivalents for the derivatives from the appendix in 
equation (28), one obtains the following equation for any arbitrary station i: 

* F,yi.2 * Giy,+3 = 1-1. 2. •••." dD 


where 
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A, = 


B.= r - 

^ h^ 12h^ ’ ^ 8h^ 6h^ 12h 


8h^ 6h 


13b. 39a. 16c. 

' ^ + 1 • 

2 ’ 




D, = 


56a. 30c, 4Cj 39a, 13b, 2a, b, 


1 


a. b. 

G,- V ■'■ J H. = 12pL (1 - eih)(l - 6ih)/Et^ 

6h^ 8h'^ ^ 0 


r(32) 


a.. = (1 - pih)(l - 6ih)' 


b^. = -2 b( 1 - 6ih)^ - 66(1 - Bih)(l - 6ih)^ 


c^. = 6 b 6(1 - 6ih)^ + 66^(1 - Bih)(l - 6ih) 


J 


In the case of first order central differences, the fictitious stations 

y-l» yn+1. and yp+p can directly be eliminated by using the finite difference 
equivalents of the Boundary conditions given by equations (30). 

In the present case of second order central differences, additional ficti- 
tious stations y_p and yp +3 must be eliminated in addition to those men- 
tioned earlier. This is accomplished by again using the symmetry conditions 

y-1 = yp y-2 = y2 (33) 


1 - fictitious stations yn+l> yp+2» yn+3 accom- 

plished by following the conditions obtained from the first order theory in the 
following manner: 


K = 0 to = -i'n-l ^ 

y“ . 0 leads to - 4y„_, + 4y„ 


(34) 

(35) 


These equations can be written in the following alternative forms: 
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(36) 


^V2 -^n-2^ - Vp (^7) 

and by recursion, one can eliminate yp +3 from 

(yn+3 “ -yp-S^ “ ®^'^n “ ■^n-1^ 

If yn+l> yn+2» •yn’*'3 f'"°"’ equations (36) to (38) are introduced into the 

second order central difference equivalent of yp" given in equation (A 8 ), one 
can show that y^' = 0. By using equations (36) and (37) in the finite differ- 
ence equivalent of y^ and setting the result to zero, one obtains the follow- 
ing relation; 

^n-2 - 2Vl * ^ ° 

Equation (39) states that the deflection at the (n - l)f^ station is the av- 
erage of the deflections at the preceding and succeeding stations. Thus, the 
deflection curve near the tip of the cantilever beam assumes a straight line 
form. Since the bending moment at the free end is zero, the condition of con- 
stant slope near the tip is justified and, thus, the boundary conditions repre- 
sented by equations (36) to (38) should give accurate results for n suitably 
large. 

The eigenvalue problem that results by using the second order central dif- 
ferences in equation (31) together with equations (36) to (38) is as follows: 

[R]{y) - («) 

where 
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0 


[S] 


0 


(41b) 


and 


{y} = [yiy2 • • • yp} 

RESULTS AND DISCUSSION 

The eigenvalue problems given by equations (8), (16), and (40) for the 
representative cases were solved on an IBM 370 computer by using the IMSL rou- 
tine EIGZF. The theoretical results obtained thus with the first and second 
order central difference theories are compared with the results available in 
the literature and these are presented in what follows. 


Axial Vibrations 

The following numerical data were used to determine the axial vibratory 
modes of cantilever beams having a length of 0.1524 m (6 in.), various taper 
parameters, a thickness ratio tQ/bQ = 0.20, and an aspect ratio L/Uq = 2.4. 
Table I shows the relative convergence rates of the first and second order cen- 
tral difference theories for the case of a uniform beam and one tapered beam 
with 6 (or 6) = 0.6 and a (or p) = 0.8. It can be seen from this table 
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that with n = 5, the second order central difference theory produces the fun- 
damental frequency of the uniform beam to within 0.005 percent error and the 
second mode frequency to within 0.4 percent error. A solution with at least 
n = 15 is required in the case of the first order central differences to at- 
tain this accuracy. The convergence of higher modes is, however, not so rapid 
in the second order theory with n = 5. Comparing the second order central 
difference solution with n = 10 and the first order solution with n = 25, 
one can observe that the maximum error encountered in the fourth mode frequency 
is of the order of 0.7 percent using second order theory with n = 10, while 
the corresponding error with n = 25 and using first order theory is about 0.8 
percent. From these results, it can be seen that the convergence is very rapid 
in the case of second order central differences and that the lowest five axial 
mode frequencies can be obtained to within 0.027 percent error with a beam di- 
vided into 30 segments (n = 30). The first order theory produces results with 
errors of the order of 1.0 percent, and Richardson's extrapolation of first 
order theory with n = 10, 20, 30 produces accurate results. 

Similar trends of convergence are observed for the case of tapered beams 
also, as can be seen from table I. It may be noted here that the percentage 
error magnitudes are relatively higher than the corresponding values obtained 
for the uniform beam. Accurate results up to the fifth mode can be obtained 
for all the taper parameters studied here by using a second order theory with 
n = 30 for all practical purposes. The frequency ratios presented in table II 
can be considered as close lower bound solutions. A further comparison of the 
first two axial mode frequencies with those presented in reference 17 indicates 
a close agreement. 


Torsional Vibrations 

In order to study the torsional vibrational characteristics, use is made 
of the same numerical data employed in the case of axial vibrations. As is 
well known, the torsional equation of motion can be obtained by replacing E 
in the equation for axial motion by (C/Ip), and consequently, the convergence 
patterns for the torsional vibrational frequencies for the case of uniform beam 
are identical to the corresponding values presented in table I for the case of 
axial motion. Since the convergence trends have already been established by 
the axial vibratory motion study, the torsional frequency parameter ratios (de- 
fined as the ratio of the square of the natural torsional frequency of any mode 
to the square of the fundamental torsional frequency) were obtained by using 
the second order central difference theory with n = 30. The ratios for the 
lowest five torsional modes are presented in figures 1 to 5. Comparison of the 
present theoretical frequencies with results obtained from the Holzer method 
(ref. 18) and the Reissner method (ref. 20) indicates close agreement. 

Table III shows a comparison of the fundamental torsional characteristic 
function of a uniform beam obtained from the second order central difference 
theory with n = 30 and the exact characteristic function. It can be seen 
that the theoretical results agree with the exact solution (ref. 19) up to five 
significant digits. Further, the higher mode characteristic functions are also 
in extremely close agreement with the corresponding exact values. 

The numerical data given in references 18 and 21 were used to determine 
the torsional frequencies of pretwisted cantilever beams. These frequencies 
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were modified by applying a correction factor given by Carnegie (ref. 16) to 
account for the fiber bending effects, and these results are presented in fig- 
ure 6. Second order central difference theory was used with n = 15. It can 
be seen that the present results agree closely with the experimental results 
of Carnegie (ref. 16) and with the theoretical results from the Reissner 
method (ref. 20). Figure 7 shows a comparison of the present natural frequen- 
cies of pretwisted tapered blading with experimental results (ref. 22) and the 
theoretical results obtained by using the Galerkin technique. A correction 
factor, applicable to uniform beams (ref. 16), was applied to the torsional 
frequencies obtained here by using second order central difference theory with 
n = 15. Close agreement between the various methods can be seen here also. 


Flexural Vibrations 

The following numerical data were used to study the flexural vibrations 
of cantilever beams having a length of 0.254 m (10 in.) and various breadth and 
thickness taper ratios: thickness ratio tp/b^ = 1.0; aspect ratio L/b^ = 40. 

The shape factor for this uniform beam ^ I /AL^ is 0.00721688 so that the 
higher order effects like shear deformation and rotary inertia can safely be 
ignored. The breadth and thickness taper parameters were varied from -0.75 to 
0.75 in steps of 0.25, and various combinations of these taper parameters were 
studied. 

A convergence study has been made for the case of a uniform beam with the 
beam divided into an odd number of segments (n = 5, 11, 15, 17, 23, 25, 29) and 
an even number of segments (n = 6, 10, 12, 18, 20, 24, 30). Both first and 
second order central difference theories were used, and the frequency ratios 
are shown in table IV. These values are also shown in graphical form in fig- 
ure 8; it can be seen from this figure that the convergence is monotonic from 
below for both the first and second order central differences and that the con- 
vergence is continuous for even and odd values of n. 

The convergence rates of the two methods can be compared from table IV or 
figure 8. As has been observed in the earlier cases of axial or torsional vi- 
brations, the lowest two flexural mode frequencies given by the second order 
theory with n = 5 are better than those given by the first order theory with 
n = 10. With n = 30, the second order central difference theory produces the 
lowest five flexural frequencies to within 0.2 percent error while the first 
order theory shows errors of the order of 2.6 percent. Using the second order 
theory with n = 30, the frequency parameter ratios (defined as the ratio of 
the square of the natural frequency of any mode of a tapered beam to the 
square of the fundamental flexural frequency of a uniform beam of comparable 
dimensions at the root) are calculated and presented in figures 9 to 19 for 
the lowest five flexural modes. Graphs are drawn showing the effects of 
breadth (or depth) taper for a given depth (or breadth) taper. Comparisons 
are made with the theoretical results, obtained earlier by using the Reissner 
method (ref. 23) and the Galerkin process (ref. 24), and also with experimen- 
tal results (refs. 11 and 25). Close agreement of the results obtained by the 
various approaches is observed. The characteristic functions obtained by 
using the first and second order theories are presented in table V for the case 
of a uniform beam where further comparisons are made with the exact character- 
istic functions (ref. 19). It has been observed that the second order central 
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difference theory gives characteristic functions close to the exact ones. 
Figures 20 to 24 show the mode shapes of tapered cantilever beams obtained by 
means of the second order theory. These mode shapes agree very closely with 
the experimental mode shapes obtained by Carnegie and Dawson (ref. 25). 


CONCLUDING REMARKS 

The second order finite difference method has been successfully applied 
to determine the uncoupled dynamic characteristics of cantilevered beams hav- 
ing variable mass and elasticity properties. Simple recursive relations have 
been used to eliminate the fictitious stations outside the beam domain by mak- 
ing logical extensions from the first order theory. The present approach is 
shown to produce accurate natural frequencies and mode shapes. The present 
improved finite difference method has the following specific advantages com- 
pared to the classical approach of using the first order central differences: 

1. For the same mesh size (step length h), the second order finite dif- 
ference method produces natural frequencies with greater accuracy than the 
first order theory. 

2. The convergence of the lower mode frequencies is very rapid in the case 
of second order central differences compared to that of the first order theory. 

3. Because of the rapid convergence shown by the present approach, accu- 
rate natural frequencies and mode shapes can be obtained directly by using a 
suitable number of segments without any necessity of the extrapolations that 
are customary with the first order central difference theory. 

4. Finally, there are few methods which produce close lower bounds; the 
present technique may be invaluable in obtaining close lower bound solutions 
without requiring extrapolations. It may be noted that even though the finite 
difference method produces lower bound solutions in general, an extrapolated 
result obtained by using the Richardson method does not necessarily give a 
bound. Thus, the present improved approach eliminates most of the shortcomings 
associated with the conventional approaches. 

The method developed in this report has the potential for extension to complex 
blade vibration problems involving coupling between in-plane and out-of-plane 
bending and torsional motions. Further extension to plate theory may prove 
beneficial since rapid convergence in the two-dimensional case may reduce the 
computational space and time considerably. 
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APPENDIX - FINITE-DIFFERENCE EQUATIONS FOR DERIVATIVES 
First Order Central Differences: 

<p'-j = ^'f’i-i “ ''' ‘i’i+1^ 

h 

q,'." = i (-q >._2 + 2<p + (P .^. 2 ) (A3) 

Zh 

(py = ^ (9i_2 - 49j_i ■*■ + (p^.^ 2 ) 

h 

Second Order Central Differences: 

‘f’i ^‘‘’i-2 “ ^i-1 ^ ^i+1 “ (A5) 

(pi = (“<p-j_2 ■'■ I®*?!-! “ “ <P-j+2)/^^*^^ (^^^ 

cp'l = (<P-j_3 “ ®‘P-j_2 ®‘(’i+2 ~ ‘P-j+3)/^^ (A7) 


cpi.v = (-^._3 + 12,p._2 - 39q,._j + 56cp. - 39cp.^^ + 12cp.^,2 - <Pi+3)/6h'* (A8) 
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TABLE I. - RELATIVE CONVERGENCE RATES OF FREQUENCY RATIOS* OF AXIAL VIBRATION FOR TYPICAL TAPER PARAMETERS 


n 

First mode 

Second mode 

Third 

mode 

Fourth mode 

Fifth 

mode 


First 

order 

Second 

order 

First 

order 

Second 

order 

First 

order 

Second 

order 

First 

order 

Second 

order 

First 

order 

Second 

order 




Breadth taper 15 = 

0; depth taper 6 = 

0 




5 

10 

15 

20 

25 

30 

0.99589 

.99897 

.99954 

.99974 

.99984 

.99989 

0.99994644 
.99999702 
.99999941 
.99999987 
.99999996 
1 .00000000 

2.89019 
2.97232 
2.98768 
2.99307 
- 2.99556 
2.99692 

2.98783 

2.99920 

2.99984 

2.99995 

2.99998 

2.99999 

4.50158 

4.87248 

4.93408 

4.96794 

4.97946 

4.95873 

4.86226 

4.98999 

4.99796 

4.99935 

4.99973 

4.99987 

5.67232 

6.65266 

6.84432 

6.91218 

6.94372 

6.96089 

6.37886 

6.94877 

6.98929 

6.99654 

6.99857 

6.99931 

6.28782 

8.26903 

8.67058 

8.81380 

8.88056 

8.91696 

7.23831 

8.83121 

8.96348 

8.98805 

8.99503 

8.99758 

Exact 
solution 
(ref. 19) 

1 

.0 

3.0 

5.0 

7.0 

9.0 

Br 

1 — , 

eadth taper 0 = 0.6, 

depth taper 6 = 0.8 

or breadth taper 

B = 0.8, depth taper 

6 = 0.6 


5 

10 

15 

20 

25 

30 

1.47717 

1.50305 

1.50785 

1.50955 

1.51032 

1.51074 

1.49710 

1.50851 

1.51035 

1.51097 

1.51124 

1.51139 

3.06602 

3.23515 

3.26716 

3.27840 

3.28362 

3.28645 

3.20418 

3.27744 

3.28681 

3.28966 

3.29089 

3.29154 

4.50158 

5.01479 

5.11520 

5.15069 

5.16718 

5.17615 

4.88551 

5.15540 

5.18261 

5.18970 

5.19249 

5.19387 

5.57925 

6.72215 

6.95531 

7.03833 

7.07701 

7.09808 

6.26438 

7.04261 

7.11627 

7.13296 

7.13885 

7.14153 

6.19245 

8.28857 

8.74112 

8.90366 

8.97962 

8.02107 

7.10365 

8.86965 

9.05151 

9.09013 

9.10270 

9.10801 


Natural frequency of tapered beam in any mode 
*Frequency ratio = frequency of uniform beam In fundamental mode 





TABLE II. - COMPARISON OF FREQUENCY RATIOS* FOR TAPERED BLADING IN AXIAL VIBRATIONS: 


FINITE DIFFERENCE SOLUTIONS; n = 30 


Breadth 

taper, 

13 

Depth 

taper, 

6 

First mode 

Second mode 

Third mode 

Fourth mode 


mode 

First 

order 

Second 

order 


Second 

order 

First 

order 

Second 

order 

First 

order 

Second 

order 


Second 

order 

0 

0 

0.99989 

1 .00000 

2.99692 

2.99999 

4.98573 

4.99987 

6.96089 

6.99931 

8.91696 

8.99758 

0 

.2 

1.04543 

1.04557 

3.01268 

3.01502 

4.99504 

5.00916 

6.96733 

7.00518 

8.92176 

9.00002 

0 

.4 

1.01460 

1.10481 

3.03616 

3.03950 

5.00911 

5.02368 

6.97713 

7.01614 

8.92908 

9.01044 

0 

.6 

1.18625 

1.18653^ 

3.07709 

3.08070 

5.03462 

5.04963 

6.99510 

7.03472 

8.94257 

9.02467 

0 

.8 

1.31006 

1.31047 

3.16917 

3.17333 

5.09932 

5.11536 

7.04296 

7.08405 

8.97936 

9.06332 

.2 

.2 

1.09189 

1.09208 

3.02991 

3.03322 

5.00528 

5.01980 

6.97444 

7.01340 

8.92707 

9.00835 

.2 

.4 

1.15215 

1.15240 

3.05523 

3.05872 

5.02055 

5.03537 

6.98509 

7.02445 

8.93508 

9.01680 

.2 

.6 

1.23511 

1.23545 

3.09858 

3.10236 

5.04772 

5.06303 

7.00427 

7.04430 

8.94943 

9.03204 

.2 

.8 

1.36046 

1.36093 

3.19405 

3.19843 

5.11506 

5.13147 

7.05413 

7.09573 

8.89777 

9.07238 

.4 

.4 

1.21366 

1.21397 

3.08289 

3.08659 

5.03739 

5.05256 

6.99687 

7.03700 

8.94384 

9.02620 

.4 

.6 

1.29802 

1.29843 

3.12938 

3.13340 

5.06679 

5.08251 

7.01768 

7.05827 

8.95949 

9.04280 

.4 

.8 

1.42465 

1.42521 

3.22938 

3.23404 

5.13781 

5.15474 

7.07036 

7.11269 

9.00003 

9.08554 

.6 

.6 

1.38368 

1.38418 

3.18019 

3.18459 

5.09947 

5.11584 

7.04094 

7.08243 

8.97702 

9.06146 

.6 

.8 

1.51074 

1.51139 

3.28645 

3.29154 

5.17615 

5.19387 

7.09808 

7.14153 

9.02107 

9.10801 

.8 

.8 

1.63529 

1.63605 

3.40099 

3.40678 

5.26289 

5.28214 

7.16392 

7.20971 

9.07223 

9.16222 


*Frequency ratio = 


Natural frequency of tapered beam In any mode 
Natural frequency of uniform beam in fundamental mode 












TABLE III. - COMPARISON OF CHARACTERISTIC FUNCTIONS OF AXIAL OR TORSIONAL 
VIBRATIONS OF UNIFORM CANTILEVER BEAM: SECOND ORDER 


CENTRAL DIFFERENCE THEORY; n » 30 


Axial 

fractional 
length, 
n = 2 /L 

First mode 

Second mode; 
theoretical 

Third mode; 
theoretical 

Fourth mode; 
theoretical 

Fifth mode; 
theoretical 

Exact 
solution 
(ref. 19) 

Theoretical 


0.00000 

0.00000 

-0.00000 

0.00000 

0.00000 

0.00000 


.15643 

.15643 

-.45399 

-.70711 

-.89101 

-.98769 


.30902 

.30902 

-.80902 

-1.00000 

-.80902 

-.30902 


.45399 

.45399 

-.98769 

-.70711 _ 

.15643 

.89101 


.58779 

.58779 

-.95106 

0.45157x10-13 

.95106 

.58779 


.70711 

.70711 

-.70711 

.70711 

.70711 

-.70711 


.80902 

.80902 

-.30902 

1.00000 


-.80902 


.89101 

.89101 

.15643 

.70711 


.45399 

.8 

.95106 

.95106 

.58779 

-.47591x10-1^ 

-.58779 

.95106 

.9 

.98769 

.98769 

.89101 

-.70711 ; 

.45399 

-.15643 

1.0 

1 .00000 

1 .00000 

1.00000 

-1.00000 

1 .00000 

-1.00000 











TABLE IV. - CONVERGENCE PATTERN OF FLEXURAL FREQUENCIES OF UNIFORM CANTILEVER BEAM USING FIRST 


AND SECOND ORDER CENTRAL DIFFERENCE THEORIES; NONDIMENSIONAL FREQUENCIES \>J pAL^EI 


n 

First mode 


Third mode 

Fourth mode 

Fifth mode 


First 

Second 

First 


DSI 

Second 

First 

Second 


Second 


order 

order 

order 


IhI 

order 

order 

order 

BeII 

order 

5 

3.4021 

3.5062 

18.870 

21.204 

45.334 

54.421 

72.469 

91.02 

92.603 

118.85 

6 

3.4359 

3.5104 

19.107 

21.553 

49.370 

57.320 

83.386 

101.85 

114.54 

144.60 

10 

3.4866 

3.5148 

21.134 

21.934 

56.603 

60.774 

104.52 

116.59 

160.58 

185.93 

n 

3.4916 

3.5152 

21.284 

21.960 

57.423 

61.014 

107.06 

117.71 

166.45 

189.46 

12 

3.4955 

3.5153 

21.340 

21.977 

58.063 

61.179 

109.07 

118.48 

171.13 

191.94 

15 

3.5029 

3.5158 

21.623 

22.006 

59.318 

61.443 

113.07 

119.72 

180.62 

195.99 

17 

3.5057 

3.5158 

21.713 

22.015 

59.827 

61.526 

114.72 

120.12 

184.59 

197.30 

18 

3.5069 

3.5158 

21.750 

22.018 

60.024 

61.555 

115.35 

120.25 

186.14 

197.74 

20 

3.5086 

3.5159 

21.801 

22.023 

60.334 

61.595 

116.37 

120.44 

188.62 

198.36 

23 

3.5104 

3.5159 

21.857 

22.027 

60.660 

61.632 

117.44 

120.61 

191.25 

198.92 

24 

3.5108 

3.5159 

21.872 

22.028 

60.744 

61.640 

117.72 

120.65 

191.93 

199.04 

25 

3.5113 

3.5159 

21.884 

22.029 

60.817 

61.647 

117.96 

120.68 

192.53 

199.15 

29 

3.5125 

3.5160 

21.923 

22.031 

61.041 

61.666 

118.70 

120.76 

194.37 

199.42 

30 

3.5127 

3.5160 

21.930 

22.031 

61.083 

61.669 

118.85 

120.78 

194.72 

199.47 

Exact 
solution 
(ref. 19) i 

3.5160 

22.0345 

61.69073 

120.9010 

199.8604 

Percent ’ 
error 
based on 
n = 30 : 

-0.094 

0.000 

-0.474 

-0.016 

-0.996 

-0.046 

-1.696 



-0.195 
































TABLE V. - COMPARISON OF CHARACTERISTIC FUNCTIONS OF FLEXURAL 


VIBRATION OF UNIFORM CANTILEVER BEAM: n = 30 


Axial Tract lorial 

Exact 

First order 

Second order 

length. 

solution 

central 

central 

n = 2 /L 

(ref. 19) 

differences 

differences 

First mode 

0.0 

0.000000 

0.0 

0.0 

.1 

.016775 

.016848 

.016783 

.2 

.063870 

.063992 

.063887 

.3 

.136485 

.136629 

.136508 

.4 

.229885 

.230038 

.229919 

.5 

.339525 

.339671 

.339565 

.6 

.461135 

.461267 

.461181 

.7 

.590875 

.590985 

.590922 

.8 

.725480 

.725557 

.725517 

.9 

.86240 

.862443 

.862424 

1.0 

1 .000000 

1 .000000 

1 .000000 

Second mode 

0.0 

0.0 

0.0 

0.0 

.1 

.092630 

.094109 

.092979 

.2 

.301056 

.303290 

.301818 

.3 

.526135 

.528641 

.527350 

.4 

.683470 

.685965 

.685029 

.5 

.713665 

.715997 

.715351 

.6 

.589475 

.591552 

.591031 

.7 

.317050 

.318772 

.318252 

.8 

-.070035 

-.068810 

-.069318 

.9 

-.523750 

-.523154 

-.523497 

1.0 

-1.000000 

-1.000000 

-1.000000 

Third mode 

0.0 

0.0 

0.0 

0.0 

.1 

.228070 

.234307 

.230020 

.2 

.604505 

.612659 

.608300 

.3 

.756240 

.763252 

.760600 

.4 

.525925 

.530052 

.528838 

.5 

.019685 

.020298 

.019714 

.6 

-.473765 

-.476490 

-.476608 

.7 

-.657425 

-.662509 

-.661663 

.8 

-.394875 

-.400383 

-.398434 

.9 

.228510 

.225100 

.226988 

1.0 

1 .000000 

1 .000000 

1 .000000 


Fourth 

mode 


0.0 

0.0 

0.0 

0.0 

.1 

.385010 

.400286 

.390900 

.2 


.768408 

.762570 

.3 

.433870 

.437896 

.438329 

.4 

-.315560 

-.322786 

-.319340 

.5 

-.707120 

-.719090 

-.715030 

.6 

-.326495 

-.334367 

-.330318 

.7 

.397390 

.399734 

.401725 

.8 

.643040 

.654228 

.651489 

.9 

.052035 

.062183 

.057070 

1.0 

1 .000000 

1 .000000 

1 .000000 

Fifth mode 

0,0 

0.0 

0.0 

0.0 

.1 

.537245 

.565373 

.550025 

.2 

.659625 

.673812 

.671266 

.3 

-.211285 

-.224220 

-.216402 

.4 

-.696550 
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Figure 8. - Relative convergence pattern of fundamental flex 
ural frequency. 





Frequency parameter ratio 



Figure 10. - Frequency parameter ratios of tapered canti- 
levered beams vibrating in flexure. First mode - effect 
of depth taper. 
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Figure 17, - Frequency parameter ratios of tapered canti- Figure 18, - Frequency parameter ratios of tapered canti- 
levered beams vibrating in flexure. Fourth mode - effect levered beams vibrating In flexure. Fifth mode - effect 

of depth taper. of breadth taper. 
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Figure 19. - Frequency parameter ratios of tapered cantilevered 
beams vibrating in flexure. Fifth mode - effect of depth 
taper. 



Figure 20. - Mode shapes of tapered cantilevered beams vi- 
brating in flexure. First mode. Breadth taper p = 0. 
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A Exact solution (ref. 19) 



Figure 23. - Mode shapes of tapered cantilevered beams vi- 
brating in flexure. Fourth mode. Breadth taper p = 0. 



Figure 24. - Mode shapes of tapered cantilevered beams vi- 
brating in flexure. Fifth mode. Breadth taper p= 0. 
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